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^p H ' Abstract 

(—i | Hard scattering processes involving hadrons at small x are described by a fey-factor- 

ization formula driven by a BFKL gluon. We explore the equivalence of this description 
to a collinear-factorization approach in which the anomalous dimensions ^ gg and 'jqg/as 
are expressed as power series in as log(l/x), or to be precise as/u> where u> is the moment 
index. In particular we confront the collinear-factorization expansion with that extracted 
from the BFKL approach with running coupling included. 
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Recently there have been several studies [1-5] of the validity and possible modification of 
the conventional Altarelli-Parisi (or GLAP) description of deep inelastic scattering in the small 
x region that has become accessible at HERA, x ~ 10 -4 . The relevant modifications are the 
inclusion of contributions which are enhanced by powers of log(l/x), but which lie outside the 
leading (and next-to-leading) Altarelli-Parisi perturbative expansion. Formally they correspond 
to the expansion of the anomalous dimensions "f gg and ■jqg/as as power series in a s /w where 
uj is the moment index. An alternative approach which automatically resums all these leading 
log(l/x) contributions to 7 99 and •jqg/as is provided by the BFKL equation coupled with the 
/^-factorization formula for calculating observable quantities || [jj. The main aim of this 
paper is to explore the connection between these two approaches. To be specific we study the 
relation between the collinear-factorization formula with log(l/a;) terms included and the kx- 
factorization formula based on the solution of the BFKL equation [Q] with running coupling as- 
We show that both approaches generate the same first few terms in the perturbative expansion 
of 7 9g and, more important, of 7 gs , which are presumably the most relevant contributions for the 
description of deep inelastic scattering in the HERA range. They differ substantially, however, 
in the asymptotically small x regime. 

Deep inelastic unpolarised electron-proton scattering may be described in terms of two 
structure functions, F 2 (x, Q 2 ) and Fl(x,Q 2 ). As usual, the kinematic variables are defined 
to be Q 2 = —q 2 and x = Q 2 /2p.q, where p and q are the four-momenta of the incoming 
proton and virtual photon probe respectively. At small values of x, x < 10 -3 , these observables 
reflect the distribution of gluons in the proton, which are by far the dominant partons in this 
kinematic region. The precise connection between the small x structure functions and the gluon 
distribution is given by the A; T -factorization formula || |7], 

F^Q 2 ) = J d -^ £^^(^4,Q 2 )/(x',4) (1) 

with % = 2, L, which is displayed pictorially in Fig. 1. The gluon distribution f(x, k 2 -,), uninte- 
grated over k\, is a solution of the BFKL equation, while F^ 9 are the off-shell gluon structure 
functions which at lowest-order are determined by the quark box (and crossed-box) contribu- 
tions to photon-gluon fusion, see Fig. 1. 

For sufficiently large values of Q 2 the leading- twist contribution is dominant, and it is 
most transparent to discuss the Q 2 evolution of Fi(x,Q 2 ) in terms of moments. Then the x' 
convolution of ([!]) factorizes to give 

^ Q 2 ) = I k 2 T , Q 2 )!^, k 2 T ) (2) 

where the moment function 

/(a,, 4)= ( 3 ) 

JO X 

with similar relations for Fi and Fj g . 
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Fixed as : /^-factorization to collinear- factorization 

It is illuminating to first consider the case of fixed coupling as [fj]. Then the photon-gluon 
moments, F~J 9 are simply functions of r = Q 2 jk\ (and uS), for massless quarks. Hence (fj) 
becomes a convolution in k\ which, in analogy with the x' convolution, may be factorized 
by taking moments a second time. In this way we obtain representations for the Fi with 
factorizable integrands 

F t (uj,Q 2 ) = — / d 7 7) /(<<>, 7)(Q 2 ) 7 (4) 

Z7TZ J c—ioo 

with c = ~. The (double) moments F? 9 and / of the gluon structure functions and the gluon 
distribution are respectively defined by 

#>, 7 ) = fdTT^TTfaT) (5) 

/>,7) = / "d4(4r 7-1 7(^,4) (6) 



where the .F 79 are dimensionless, but / carries the dimension (/eq) -7-1 . Representation ([|) en- 
ables the leading-twist contribution to be identified from a knowledge of the analytic properties 
of / and F? 9 in the complex 7 plane. 

The gluon distribution f(x, k\) satisfies the BFKL equation, which in moment space has 
the form 

/(u,, 4) = f>, 4) + ^ / ff^(4, #)7(«, 4 2 ) (7) 

where a$ = 30:5/71 and F is the usual BFKL kernel. The double-moment function / is therefore 
given by 

/>, 7 )= f ° (un l (8) 
1 - fcsMKi-y) 

where Ki^y) is the eigenvalue of the BFKL kernel corresponding to the eigenfunction propor- 
tional to (4) 7 - it can be shown that 



k{i) = 2^(1) - ^(1 -7) - ^(7) 
1 
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l + £2C(2n+l) 7 2n+l 

n=l 



(9) 



where \I/ is the logarithmic derivative of the Euler gamma function, ^f(z) = T'(z)/T(z), and 
C(n) is the Riemann zeta function. 

We see from (|]) that the large Q 2 behaviour of Fi(u, Q 2 ) is controlled by the pole at 7 = 7 
of f(u, 7) of (^|) which lies to the left of, and nearest to, the contour of integration in the 
7-plane. For a physically reasonable choice of input f°, this pole arises from the zero of the 
denominator of (§|). That is 

1 7.R 



\-{a s /u)K{i) 7-7' 



(10) 
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where from (^j) we have 

7 = 
R = 



^ + 2C(3)(^) 4 + 2 C (5)(^) 6 + o(^) 7 , (11) 

CO \ OJ J \0J J \ UJ J 




l + 6C(3)f^)Vo(^) 5 . (12) 



7 is the leading-twist anomalous dimension . If we insert the pole of ([T0|) and fl8|) into (f|) , and 
we close the contour of integration in the left-half plane, then we obtain the high Q 2 behaviour 

F r (u,Q 2 ) = F? 9 (un) Ti? (jfj f°(un)(Q 2 T- (13) 

Eq. (p~3|) is the usual formula for the factorization of collinear (or mass) singularities written in 
moment space. This becomes more apparent if we express ([P3|) in the form 

F t (co,Q 2 ) = Cj 9 (u;^)g(u J ,Q 2 ) (14) 

where _ 

^7 9 (^7) = 7^ 79 (^7) J R(^f) (15) 

is the moment of the (process dependent) coefficient function and 

g{u>,Q 2 ) = {QlVf\u ll ){^\ 1 (16) 



is the moment function of the (integrated) gluon density. Thus we can identify (Qq) 7 Pi^il) 
with the moment of the gluon distribution at the "starting" scale of the evolution in Q 2 . 

The quantity R, the residue in ([H]), is renormalisation scheme dependent [0. For studies of 
the BFKL equation, ([?[), it is appropriate to regularise R by choosing an inhomogeneous term 
of the form 

f°(oj,k 2 )=G°(oj)8(k 2 T - f i 2 ). (17) 

On the other hand 7, and F? 9 , are scheme independent (at least in the so-called regular 
schemes). The scheme dependence of R is compensated by subleading contributions of 
0(as(<ys/ UJ ) n ) i n the anomalous dimension 7 9g , which are still at present unknown. This 
cancellation takes place when we allow the coupling to run. 

The above collinear factorization formula ([H|) is true as it stands for F^, but some care is 
needed for F 2 . First we check its validity for F L . Since — ► constant for large r, we see from 
(|5|) that F1 9 {oj,^) ~ I/7. However, this potential singularity at 7 = is cancelled by the 7 
factor in the numerator of ([3~0|) . On the other hand F 2 9 log(r) and hence i*g (a>,7) ~ I/7 2 , 



3 



where the double pole reflects the collinear singularity associated with the g — > qq transition. 
The integration contour in (H) therefore also encloses the pole 

j^( W)7 )/( w , 7 )~ Jj7~-, (18) 
7 7 

which gives rise to a "scaling sea" contribution to F% which is independent of Q 2 . To remove this 
contribution and to focus attention on the effects of the perturbative pole at 7 = 7 we consider 
the observable dF 2 / d log Q 2 , rather than F 2 itself. In this case the collinear factorization formula 
is of the form 

^|^ = E2e?P^,7)^Q 2 ) (19) 
where the coefficient, or gluon-quark splitting, function is given by 

P qg (u,j) = a s ^u;,j) R ^j , (20) 

with 

a 5 $(^,7) = 7 2 i^,7) (21) 

defined to be a regular function at 7 = 0. Since $ is known [Q] in terms of the quark box (and 
crossed box), we can determine the perturbative expansion of P qg by calculating 



P qg (^,l) = as 



<9$ 



+ 

7=0 



R (22) 



with 7 and R given by the perturbative expansions of ([□]) and ( |P2| ) respectively. To be precise 
we substitute for 7 in fl2"2] ) and obtain the perturbative expansion of P qg /as as a power series 
in cts/uj. This provides the recipe to compute the leading log(l/x) contribution to P qg . 

Although we have expanded the observables in a perturbative series in as/u>, we should 
recall that the small x behaviour of Fl and dF2/d\ogQ 2 is controlled by the singularities of 
Fi(u, Q 2 ) in the u complex plane. The singularities arise from 7. The leading singularity of 7 
is the BFKL branch point at ui = uo^ = (41og2)a 7 s'. To see this we note that the position of 
the singularity is controlled by the value of K{j) at its symmetry point, 7 = |. We expand K 
about this point 

£(7) = 4 log 2 + 14C(3)( 7 - \f + . . . , (23) 
and determine the leading singularity as the implicit solution of 

1 - ^(7) = 0, (24) 



see (||). The leading pole of /(a;, 7), which lies inside the contour of integration of (f|), is at 



where u>l = (4 log 2) a s . In x space the leading singularity of the anomalous dimension gives 
an x~^ h behaviour of the gluon distribution at asymptotically small values of x. On the other 
hand the perturbation series in as in (|TTD, enables the collinear factorization formula to 

be used to investigate the approach to the BFKL x~ UJL form, as x decreases since 

V- ( Us Y \- (a s logl/x) n - 1 

The collinear factorization formulae give well-defined perturbative expansions for Fl and 
dF%jd log Q 2 which allow the leading as log(l/x) contributions to be resummed. We will 
discuss the implications of the reduction of fey-factorization to collinear form after we have 
implemented the running of as- 

Running as collinear-factorization to /^-factorization 

To see the effect of the running of as we simply replace 

j(a s ,co) -> 7(«5(Q 2 ),^), (27) 

and similarly for Riois/uj), in ([14]) and (|i~9"l). The crucial change is in the Q 2 evolution factor 
of g(to, Q 2 ), which becomes 

{%S - AC ■ (28) 

In the small x BFKL limit 7 is simply a function of the ratio o7,s(g 2 )/u;, as in the fixed cou- 
pling case. We see immediately the important role played by the non-perturbative region. To 
illustrate the effect, it is sufficient to take 

^ f„2\ _ l/i„„/„2 / A 2 



a s {q I )=b/\og{ q 2 /k 2 ), (29) 



and to write (P~T| ) in the form 



where the coefficients are known (and in particular A\ = 1 and A 2 = A 3 = A 5 = 0). Then the 
exponent in (|28|) is given by 

Q 2 dq*_ b ( a s (Ql) \ b " A n j (asm V"' ,, n 

Q} Y^^ oe {^m) + ^—i\{—) -{—) ]■ (31) 

The first term on the right-hand-side leads to the usual double-leading-logarithmic (DLL) 
behaviour of the gluon distribution g(u, Q 2 ). The sum in the second term builds up the BFKL 
behaviour, and here we see the dominance of the asiQo)/ 00 P ower series evaluated at the 
starting scale Qq as compared to the truly perturbative power series in o~s(Q 2 )/uj. In other 
words the leading singularity in fl3i"D is the BFKL branch point at uj — a"s((5o)41og2. In 



principle it should be reabsorbed in the starting distribution g(u, Q%), leaving the perturbative 
contribution which is controlled by as{Q 2 )\ yet in practice it is the full formula ([H]) which is 
used § i g |. 

The above representation contains therefore an equivalent infrared sensitivity to that con- 
tained in the direct BFKL predictions J/J, but we see that it has been explicitly isolated in a 
factorizable form. By infrared sensitivity we mean that the leading singularity in the uj plane is 
controlled by Qs(Qo) an d not by as(Q 2 ). The distinction between Q 2 and Qo i s 5 °f course, im- 
material in the region (Q 2 > Qo) °f applicability of the genuine leading log(l/x) approximation, 
that is a s {Ql) \og{Q 2 /Q 2 ) « 1, but a s (Q 2 ) \og(l/x) ~ O(l). 

The BFKL corrections to the DLL contribution only enter the expansion for the anomalous 
dimension (|3TD at order (as/uj) 4 and above, whereas for P qg of (^) it can be shown that all 
terms (n = 0,1,...) are present in the expansion B n (as/uj) n . For this reason we expect 
that the small x behaviour of dF 2 /dlogQ 2 in the HERA regime will be controlled more by the 
perturbative expansion of P qg than of 7. However, as x decreases the expansion of 7 will begin 
to play a dominant role. 

The BFKL equation was originally derived for fixed as- The correct way to include the 
running of as is not firmly established. The procedure usually adopted is to take as{k^) 
in (|7|) so that the DLL limit of GLAP evolution is obtained. Here we find the perturbative 
expansion obtained from this prescription. We are therefore able to check the validity of 
the procedure by comparing with the expansion obtained from the renormalization group (or 
collinear factorization) approach, that is (|27| ) -(|3T| ) . 

If we replace the fixed as of (0) by as(k^) the BFKL equation becomes 

7(", 4) = log f>, k 2 T ) + h - J ^K(k 2 T , k> 2 )J(u, 4), (32) 

which, in terms of the moment variable 7 conjugate^ to k^/A 2 , reduces to the differential 
equation [10-13] 

cry cry uj 

From the extension of @) and (f21~l) to running as we see that the fc^-factorization gives 

dF 2 {u,Q 2 ) 1 ri-Hoo _ 1 /gay 

and similarly for Fl(u,Q 2 ), where the double-moment of the gluon f{oJ, 7) is the solution of 
( |33| ) . The leading-twist contribution is controlled by the solution of the homogeneous form of 




(II) [11, 12 



b 

UJ J-y 



dik{i) 



(35) 



2 For running as we choose to inter-relate dimensionless quantities /<->/, whereas for fixed as it was 
convenient to allow / to carry dimensions, see (^|). 
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where H°(u) will eventually have to be fixed by the starting gluon distribution, f°. 

We now inspect the perturbative expansion of ([34] ) and find that the first few terms are 
identical to those in the renormalization group expansion. To make this identification we 
concentrate on the expansion in terms of powers of as{Q 2 )/^ which contain the hard scale Q 2 . 
The non-perturbative contributions can always be absorbed into a redefinition of H°(uj). 

To begin we note that the leading twist contribution is controlled by the strip — 1 < 7 < 
of the branch cut in the 7 plane, where the branch point at 7 = is generated entirely by the 
I/7 term in K(i). We isolate this singularity by introducing the function 

F( 7 ) = K(i) - - (36) 
7 

such that K is regular at 7 = 0. We then insert the result 



- / di Kirf) = log 7 + - / t/y K{i) (37) 

into (35) and (33), fold the contour around the cut and evaluate the discontinuity to obtain 

9F }} U ^P = ~ sin (—) H\u)a s {Q 2 )I + higher twist (38) 
a log Q 2 \ uj J 

where ^ 

j = |° i rf 7 $(c,7)(-7)-^ 1 exp r- J^diK{i) \ (^A . (39) 

The integral I has, of course, to be understood in the sense of an analytic continuation since 
it diverges at 7 = 0. To expand J in a perturbation series we first change the variable of 
integration 

7 -> -p/t where t = log(Q 2 / A2 ) ; (40) 
then the integral takes the form 

rt . , (b_ 

^UJ J-r/i 

We use (0) to expand the first exponential factor in (f41|) 



t b/uJ f dp $(lu, -p/t)p-^- 1 exp [ - / drf K(i)) e~ p . (41) 
Jo \u J-p/t I 



where we have omitted a factorizable non-perturbative contribution coming from the upper 
limit. When we expand the exponential in ([E3) and insert the series into ( flTD we encounter 
integrals of the form 



dpp-^ +2n e- p = r(-- + 2n + l 



\2n+l. 



r(-6/w) (-b/uY n+l (l + 0{u)). (43) 
7 



Here the contribution from t < p < oo gives higher-twist terms which vanish as 1/Q 2 , modulo 
logarithmic corrections. The term T(—b/oj) can be reabsorbed into the starting distribution, 
where it belongs, and we find the perturbative expansion of / is of the form 



t b /"T(-b/uj) $(w,0) 



1 - 



oj 



E 

n=l 



2C(2n+l) (a s {Q 2 



2n + 1 



OJ 



+ higher order terms. (44) 



We see that the first two terms (n = 1, 2) are identical to the first two terms (n = 4, 6) in the 
perturbative expansion of ((£[]), which are proportional to as(Q 2 ) 3 an d as(Q 2 ) 5 respectively. 
The DLL contribution in fl3T|) corresponds to the t h ^ factor in (H). 

To generate the expansion of P qg /ots as a power series in as/oj we expand the function 
$(u;,7) of (JT) around 7 = 0. This procedure generates the same first three terms as those in 
the expansion shown in (|22|). At higher order, (a~s/oj) 3 and above, we see that the terms of 



O(oj) in (f43|), as well as various other contributions, will also contribute to the expansion of 
P qg . However, it is the first few terms that are important for the onset of the BFKL behaviour 
in the HERA small x regime [I], fj. Note that the perturbative expansion in ( |4"4| ) contains 
an additional factor of b/u which enables this series to be separated from the perturbative 
expansion of <3>(u;,7). 

Before we conclude, we can gain further insight into the relation between the BFKL equation 
and collinear factorization in the case of running as if we estimate the integral fl34]) using the 
saddle-point method. Inserting (|35| ) we see that the position of the saddle-point, 7, is given by 
the implicit equation 

b ^ , , (Q 2 ^ 



0. 



that is by 



as(Q 2 ) 



which is the same as (|24|) for fixed as- We evaluate the integrand of (|34j ) at 7 
to rearrange the product /(oj,^) (Q 2 /A 2 ) 7 in the form 



(45) 

7 and use (|45|) 



H\oj) exp j (± J_ dj> Ktf)\ + 7 log 



H°(oj)exp 



C ¥ T(«5(? 2 )/o;)} (46) 



where H° includes the integration constant. The equality ([46]) is obtained by integrating the 
integral on the left- hand- side by parts. Thus the saddle-point estimate of (|34"D gives 



OF 2 (oj,Q 2 



as{Q ) , _ = exp 




l(a s (q 2 )/oj) 



where 



dlogQ 2 \]-i 2 K' 

K' = dK/d^j _ and 7 = ^{ois{Q 2 ) / oj). 



(47) 



7=7 
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This representation is applicable in the region uj > ujl = (4log2)as{Qo)- F° r smaller values 
of uj the saddle-point estimate involves two (stationary phase) contributions which lead to a 
different representation of the integral. In other words ( fT7j) is not a valid approximation of the 
integral fl3~4]) for < Unlike the case of fixed as, the BFKL solution for running as does 
not contain the branch point singularity at uj = ujl, but rather it has (an infinite number of) 
poles in the uj plane. The poles are determined by the starting point condition and hence are 
controlled by as at Qq. It turns out that the leading pole singularity occurs at uj < Ul{ols{Q^)) 



14, 13 



In summary, we have confronted the collinear-factorization approach for the calculation of 
observable quantities at small x with the evaluation based on the fc^-factorization formula. 
For fixed as both approaches are equivalent at the leading- twist level. In fact the insertion 
of the solution of the BFKL equation into the /^-factorization formula provides a recipe for 
calculating 7 = 7 99 and P qg / as as power series in as/to, where to is the moment index. The effect 
of introducing a running as in the collinear-factorization formalism is summarized by fl27|) and 
(|28|) . We noted that the leading singularity in the uj plane is a branch point at uj = ujl(Qq) which 
is controlled by asiQl) rather than as(Q 2 ), c.f. d31|). That is the truly perturbative behaviour is 
hidden behind a non-perturbative contribution. In principle, the latter could be factored off and 
absorbed into the starting distribution. We then examined fc^-factorization with BFKL input 
with running as and compared the predictions with those obtained from collinear factorization 
with running as- We found the remarkable result that both the factorization prescriptions 
generate a perturbative expansion as a power series in as{Q 2 )/^ with exactly the same first 
few non-trivial terms, on top of the same DLL contribution. These terms are the most important 
perturbative contributions for the onset of the leading log(l/x) behaviour in the HERA regime. 
In practice, in both the collinear- and fc^-factorization approaches, the leading singularity in the 
u;-plane, which controls the small x behaviour, depends on Qq. The location of the singularity 
is different, however. In the first case the power series in as(Ql)/oj builds up a branch point at 
uj = ujl(Qq), whereas in the second case we generate a leading pole at a considerably smaller 
value of uj. We conclude that the truly perturbative contributions in the two approaches are 
remarkably similar, but in practice they are partially hidden by non-perturbative terms. 
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Figure Caption 

Fig. 1 Pictorial representation of the fc^-factorization formula of (|]). 
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